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1 Prelude 



At this stage, integrable one-dimensional quantum spin chains are a widely studied 
subject in theoretical physics. In fact, they are interesting models in themselves as in- 
tegrable (qualitative) description of interactions between microscopic magnets in solid 
state physics phenomena and, besides, they may be thought of as (lattice) regulari- 
sations of two dimensional quantum field theories (cf. for instance [2j, [2j): in this 
respect, they opened the road of the application of techniques, like the Algebraic Bethe 
Ansatz, which would be otherwise considerably difficult if applied directly to the in- 
finite degrees of freedom case. Finally, recent and surprising findings (see jl] and the 
subsequent development) showed that hamiltonians of some integrable spin chains co- 
incide with the mixing matrices of gauge invariant operators of M = 4 super Yang-Mills 
theories, opening a new way to test the AdS/CFT correspondence. 

Given this growing interest in the subject, we started by studying the most general 
one dimensional spin 1/2 chain, the so-called 1/2-XYZ, which describes spin 1/2 degrees 
of freedom on a lattice with anisotropic nearest-neighbour interaction. In the infinite 
volume limit and with finite lattice spacing (thermodynamic limit), a suitable scaling 
of the elliptic nome allowed the authors of Pj to obtain an integrable quantum field 
theory in the two dimensional plane, the sine-Gordon model. Therefore, it is reasonable 
to think of the 1/2-XYZ chain (in the thermodynamic limit) as an elliptic deformation 
of the sine-Gordon field theory. Moreover, in a recent paper [S] we proved that the same 
conclusion holds if we also allow the lattice spacing to go to zero in such a way to keep fi- 
nite the lattice length on a circumference (cylinder geometry). In fact, starting from the 
Algebraic Bethe Ansatz equations we wrote down the equivalent Non-Linear Integral 
Equations (NLIEs) describing the anti-ferromagnetic ground state, its first excitations, 
- i.e. spin waves propagating along the chain -, and bound states thereof, which turn 
out to be labelled by a positive integer n 0- These excitations become in the sine- 
Gordon limit the solitons/antisolitons and breathers respectively. Therefore, we may 
think of them as elliptic deformed counterparts of solitons/antisolitons and breathers, 
and definitely borrow these names from the field theory parlance. By inspecting the 
NLIE in the thermodynamic limit, we started to compute scattering factors between 
excitations [5]. Our aim was to reveal, within a clear and unique framework, the elliptic 
scattering factors which had been already proposed [HI E], without however trans- 
parent connection to physical models. In this respect, the scattering matrix between 
excitations in the repulsive regime - elliptic solitons and antisolitons - was computed 
and showed to coincide with Zamolodchikov's S-matrix 6|. The elliptic nome of such 
a matrix was found shifted with respect to the one appearing in the parameters defin- 
ing the XYZ chain, thereby proving the so-called Smirnov's conjecture. On the other 
hand, in the attractive regime the scattering factor between two first (n = 1) elliptic 
breathers was also computed [3] and showed to coincide on one side with the Lukyanov- 
Mussardo-Penati S- factor [7||H], on the other side with the elliptic structure function of 
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the Deformed Virasoro Algebra pQ. 

In this paper we will continue the calculation of S-matrices by computing all the 
possible scattering factors between two n and m elliptic breathers, with n,m > 1. 
As expected, these factors are elliptic deformations of the corresponding ones between 
sine-Gordon breathers Pj. From the mathematical point of view they may be used, 
when n — m, to define an order n Deformed Virasoro Algebra (DVA). For the well- 
known DVA (order n — 1) [I] is obtained by quantizing the Frenkel-Reshetikhin Poisson 
bracket jTOJ and all the other higher order ones {n > 1) may be thought of in the same 
manner. 



2 An outlook on the theoretical background 

The spin 1/2-XYZ model with periodic boundary conditions is a spin chain with hamil- 
tonian, written in terms of Pauli matrices <j x,y,z , 

1 N 

V = 2 D J *«+i + VX+i + J*« + x) ■ (2-1) 

71=1 

Here N is the number of lattice sites and because of the periodicity the site N + 1 is 
identified with the site 1. The three (real) coupling constants J x , J y and J z may be 
reparametrised (up to an overall multiplicative constant) in terms of elliptic functions 
with nome exp ^— Tr^j: 

J x = 1 + k sn 2 2r? , J y = 1 - k sn 2 2r? , J z = cn2r? dn2r/ , (2.2) 

where 77 is real. As for the notations used in this paper we address the reader, for 

instance, to the Section 2 of 0. As in [3], we constrain our analysis to the disordered 
regime, i.e. 

< exp f-^J < 1 > 

0<r/<K, (2.3) 

which are equivalent to k > 0, K > 0, K' > or | J z \ < J y < J x . We also distinguish 
between the repulsive regime < 77 < y, and the attractive regime y < r\ < K. 

The spin 1/2-XYZ chain enjoys many interesting limit cases. If K' — > 00, which 
implies K — > |, we have J x — ► J y — > 1, J z — > cos2?7, so that we recover the XXZ 
chain in the massless anti-ferromagnetic regime. If K — > 00 (K 7 — > | ), we discover that 
Jy —> J z cos h^ §~ ' ^ x ~^ 1+tanh 2 2i], which means that we obtain the XXZ chain in the 
anti-ferromagnetic massive regime. As mentioned in the prelude, particular attention 
should be payed to the double scaling limit of site number and coupling constant 

877 

iV -> 00 , K' = — log N + finite terms , (2.4) 

7T 
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on the cylinder of constant circumference R and lattice spacing A = f -»• 0. In this 
case we obtain the sine-Gordon model on a cylinder of circumference R with coupling 
constant (3 2 = 1 — ^2, i.e. with Lagrangian C = \(d<p) 2 + £^ (cos \/87r/30 — 1) jH]. 



87T/3 2 

Eventually, if even R goes to infinity, we gain the plane geometry, where the famous 
results of j2] apply. 

The spectrum of the spin 1/2- XYZ model was found in jTI] by means of the Algebraic 
Bethe Ansatz. In particular, the Bethe equations read 



H(iaj + r])Q(iaj + rf) 
H{ia.j — r))Q(iaj — rj) 



N M 



-4m V JL tt H(ia.j - ia k + 2r])Q(iaj - ia k + 2rj) 
— e 2K 



n 



fe= A H(iaj - ia k - 2r))Q(iaij - ia k - 2r]) ' 

(2.5) 

where j = 1, M and v is an integer. These equations ()2.5|) are valid when 

m\T) = 2m 2 K , (2.6) 

where mi and m 2 are integers, and when 2M = N, mod m\. The corresponding 
eigenvalues of the transfer matrix An (a) may be written in terms of the Bethe roots as 

^4 = e**>*H{ia + V )»G(ia + rjf f\ ™ + - + ^ + 

6(0)* v /; v /; H(iatj - ia)e(iaj - ia) 



(2.7) 



M 

J-J- W I in — i.ru • IHI irv — m- 



i if (ia — iaj)Q(ia — zct^ 



While writing (|2.5[ I2.7jl we borrowed the notations from [H]: in particular, all the 
elliptic functions have nome exp (— tt^). 



3 The Non-Linear Integral Equation 

In j^j we described the vacuum and the first excitations of the spin 1/2-XYZ chain 
upon transforming the Bethe equations into an equivalent NLIE, one for each state. In 
Algebraic Bethe Ansatz language the vacuum is given by v = and by filling the real 
interval ] — , with Bethe roots. Excitations over the vacuum are expressed by 
Nh (real) holes, N c close (complex) pairs, N w wide (complex) pairs and N sc (complex) 
self-conjugated roots, the vacuum obviously corresponding to Nh = N c = N w = N sc = 
(cf. jS] for details). In fact, any of these states may be described in a compact way by 
means of the following NLIE for the unknown function 

N h N c N w 

Z N {6) = NF(9) + J2x(0-h k )-J2x(0-c k )-J2xn(0-w k )~ 

k=l k=l k=l 

Nsc r **L 
- Yxn(6-s k ) + 2 / 2 " drjG(6-ri)Imln U + j^^] , (3.1) 

/ KK' L J 

k=l 
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where we assume for granted from [S] the definitions of the functions of the renormalised 
rapidity 6, Z N (6) = Z N (£0) (similarly for F and X ) and G{9) = £G (£9). Nonethe- 
less, we recall for the lazy reader the explicit formulae for the forcing term, the kernel 
and the scattering phase 

+°° f -inn 2n ( K ~ r >) 7r h U.\ rn-.T77r-.inh 2nKw 

svm 1 smn K ' I 1 Kjcosmrsmn K , ( o ^ 

F{<>) ^ n smh ^ + smh MI e • < 3 ' 2 > 

n= — oo K K 



s i nh MK^ 



G{6) ~ KK' £ 2sinh^^cosh^ e " 2 ^ ' ^ 



£(0) = / dx2TrG(x) = 

sinh« e" 2 ^ K 2 



= 5 ^ — , |Im0|<2— -2K. (3.4) 

n tTl 2 Sinh 2 "(K 7 ^ CQsh gn^r n v 

Besides, we required reference to the so-called second determination of x, xiii m that 
wide roots lie, by definition, outside the strip |Im0| < min{2K , 2K 2 /r/ — 2K}, where 
formula (|3.4|) rigorously holds. Actually, here we only need the second determination 
in the attractive regime y < r\ < K 

Xiiifi) = x(0) - X(6 - 2 i (K 2 /r ] - K) sgnlmfl) , (3.5) 

or better, spitting the upper from the lower band, we may equivalently introduce the 
two functions 

jfii(0) = Xii(0) when elm > 2K 2 /^-2K, e = ±1 . (3.6) 
It is convenient to re-express (|3.5|) in terms of the multiplicative independent variables 
( 2inr,9\ / 4vrK\ / 4n V \ 

x = ex v{~KKr)' q = exp \^r)' p = exp V k^J' (3 ' 7) 

which naturally lead to the definitions (\a\ < 1) 

oo 

(x; a) = J^[(l — xa s ) , Q a (x) = (x; a)(ax~ 1 ; a)(a; a) . (3.8) 

s=0 

The latter satisfies the multiplicative properties 

Q a (ax) = Qaix- 1 ) = -x- 1 6 a (z) . (3.9) 
Thanks to the definitions ()3.8|) we may write ()3.6|) explicitly 



xg(0)=iln 



Q p 2(xq 2 p 2 )Q p 2(x l q 2 p 2 ) 
Q p 2 (xg 1 ^ i p~ 1 2 i )9p2 (x^q^p^r) 



(3.10) 
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4 The bound states, i.e. the elliptic breathers 



In the attractive regime K/2<r^<Ka soliton and an antisoliton may form a bound 
state, which may be named - in analogy with the sine-Gordon theory - elliptic breather. 
It is represented as Bethe root configuration by adding to the Fermi-Dirac sea of the 
vacuum a string of complex (conjugated) roots 1 . After remembering |2] that in the 
thermodynamic limit N —>■ oo the breather B n (8) with (real) rapidity 9 can exist only 
if the integer n satisfies 

— |jtK < 7] , (4.1) 
n + 1 

we must consider separately the two cases of n even and odd: the breather B 2 k+2(9) is 
described by the string (9 real) 

9±iK±iK2h (1 - — ) , h = 0,1,..., A;, (4.2) 



V 

while the breather B 2 k+i(9) by the string 

K 2 / K\ 

9 + i — , 9±iK±iK(2h- 1) I 1 , h = l,...,k. (4.3) 

V V V/ 

As a consequence of the constraint (|4.1|) the previous strings contain wide roots only, 
i.e. they are entirely in the second analyticity region of Zn(9). An explanation of this 
fact will be given in a forthcoming paper and will justify the constraint (j4.1)l itself. 

We are now ready to study the collision two breathers of the XYZ spectrum. Besides 
the clear interest as an example of non-relativistic factorised (integrable) theory, this 
scattering acquires a peculiar role in a mathematical perspective, since we showed in a 
previous publication j^j that in the case of two low-lying (n = 1) breathers it furnishes, 
- as Zamolodchikov-Faddeev algebra -, the Deformed Virasoro Algebra (DVA) of [T] . In 
order to obtain the asymptotic states, we shall consider the thermodynamic (N — > oo) 
limit of the chain. 



5 Scattering factors between two elliptic breathers 

When the number of lattice sites N goes to infinity, the NLIE undergoes a significant 
simplification in that the convolution term can be neglected with respect to the other 
terms [SUES Ej. Hence, when iV — ► oo the NLIE (|3.1|) . specified for the configuration 
containing an elliptic breather -E>2fc'+2(#i) an d another breather -82^+2(^2), takes on the 

1 This string was the simplest one, namely a pair, for the low-lying breather of [Jj; besides the 
soliton/antisoliton and breather excitation exist and are described by the regular string values below 
only in the thermodynamic limit N — ► 00. 
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simplified form 



Z N (9) = NF(6) */? (e-9i + teK + ie'K2ti (l--Y\- 

h'=0e'=±l ^ \ V J J 

~ EE^ (9-e 2 + ieK + ieK2h(l--) 

h=0 e=±l ^ V V / . 



(5.1) 



Now, we need to compute Z N {6) on each root of the n even string 

K 
»7 



^(e', tf) = 0! + k'K + k'K2/i' 1 



, e' = ±l, /i' = 0,l, 



(5.2) 



and sum up over all the roots (i.e. over both e' and h'). From the definition of 
counting function Z N {6) we know (see discussion in Section 6 of \5f) that the sum 
of the scattering phases (i.e. the terms involving Xiii but not the forcing term) will 
equalise exactly the scattering phase i In SW+2, 2^+2(^12), #12 = &i — #2, between the 
elliptic breathers B 2 k'+ 2 (Qi) and B 2 k+2 (#2)- in other words, we can interpret on a 
physical ground S 2 k> ^+2,2^+2(^12) as the scattering amplitude. Technicaly, we shall pay 
attention to the positions of the roots (J5.2j) in the second analyticity strip of x, i.e. 
e'Im6 l i(e / , h') > 2K. 2 /r] — 2K. Therefore, we consider the second determination of the 
functions Xii(@) m f|5. If) : this entails that we must apply again ()3.5|) . now with new 
input Xn(@) i n the r.h.s.. We obtain the new functions Xn' e \^)n 
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i In 



— xq~ 



p 2 (xq 1 



c+e 
2 



P 



Q p 2 (x 1 q 1 "a p) 9„2 (xg 



or explicitly 



J7 



z In 



% In 



— xq 



x ( I T\d)n = xfr\o)n = i\n 

Hence, the scattering factor to compute reads 



6 p 2(x 1 )Q p 2(xq 2 p) 
p 2 (x~ 1 p)O p 2 (xq 2 ) 

1 Q p 2(xp)Q p 2(x- 1 q 2 ) ' 
Qp2(x)Q p 2(x- 1 q 2 p)\ 

Q p 2 (xqp) Q p 2 (x _ 1 q) 



-x 



e p 2(x 1 qp)Q p 2(xq) 



S 2 k' 
exp 



2fc'+2,2fc+2^12 
fc' fc 



(5.3) 

(5.4) 
(5.5) 
(5.6) 

(5.7) 



E E E E *& B) (*i2 + + OK + i(e'h> + eh)2K (l 

h'—r\ h—r\ c > — k <==+ ^ ^ 



/i'=0 /i=0 e'=± e=± 

and we can use the explicit expression (|5.3j) to write down 

k' k 



K 

V 



n 



Soy 



nnnn 

h'=0/i=Oe'=±e=± L 



-X 



6p 2 (ar V" 1 g 1 - / )Bp2 (xpj^V^V) 
9p2 (a;p-/- 1 g 1 +/)0p2 (ar V _1 5 1_/ p) . 



(5i 
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with the shorthand / = ^f- + eh + e'h' and where, with a little abuse of notation, we 
defined again x as 

2^^776*12 



x = exp 

This expression may be rearranged as 



KK' 



(5.9) 



k' k 



Sou 



2fc'+2,2fc+2^12j 



nnnn 



h'=0 h=0 e'=± e=± 

and, after performing the product over e, e', as 



1 i-f®p 2 (. x V l 1 ' 1^ 1 )@ P 2 ( ;r P ^ 1 qf +1 p) 
Q p 2(xp~f~ 1 q l+ f)Q p 2(xp 1 ~f qf~ l p) 



f-i i-f^p- 



(5.10) 



fc'+i 



2/c'+2,2/c+2l,<7l2y 



n n 

h'=-k' h=-k-l 



_! e p a(V" h "' l '3 fc+A, " 1 )©p!»(a!P" fc " fc '" 1 9 ,H " h ' +1 P) 



P 9 



6 p 2 (xp" 



-h-h'-lgl+h+h' 



')Q p 2(xp 



l-h-h'gh+h'-l 



(5.11) 



Now, the product on /i' gives 

(2fc'+2)(2fc+2) fe 



2fc'+2,2fc+2l (7 12y 



py 



n 



h=-k- 



Q p 2(xp k '~ h q~ k ' +h )Q p 2(xp~ k '~ h ~ 2 q k ' +h+2 p)Q p 2(xp -~^q 



e p2 (xp 1+fc '- /i g- fc,+Ai - 1 ) 

x Q p 2(xp- k '' h - 2 q 2 + k '+ h ) 

k'-h-l„k'+h+l 



P) 



Q p 2(xp~ k '- h - 1 q 1 + k '+ h )Q p 2(xp 1 + k, - h q- k '+ h - 1 p)Q p 2(xp k, - h q- k '+ h p) 



(2fc'+2)(2fc+2) 



Q p 2(xp k> k q k k ')Q p 2(xp k+k ' +2 q k k ' 2 )Q p 2(xp k k ' 2 q k + k '+ 2 p) 



Q p 2{xp- k - k '- 2 q k + k '+ 2 )Q p 2{xp k - k 'q k '- k )Q p 2(xp k + k '+ 2 q- 



k-k'-2 



, p . 

,/.- /.'„//-/,•„, 2 2 (xp k '~ h q h ~ k ')9 2 2 (xp~ h ~ k '~ 2 q h+k ' +2 p) 



Q p 2(xp k k 'q k ' k p) 

' e p 2(xp k '- k q k - k 'p) Q2 2 ( xp ~h-k>-2qh + k>+2}Q2 2 ( xp -h+k> q h~ki p } ■ 



n 



p) 



(5.12) 



Eventually, we can reach the final form after performing algebraic manipulations in 
order to highlight the exchange symmetry x — » x" 1 and with the positions 2k! + 2 = 
n, 2k + 2 = m 



S. 



n,rn\ u 12 i 



X 



. , , _ -I n — m m — n » „ , -1 ■ m — n n — m , . . , -1 m + n m-\-n < 

2 <a p 2{x q 2 p 2 )<~} p 2{xp ^ 2 q 2 )<S p 2{x q 2 p 2 ) 



-~ 1 II — lit lit — IL \ , -, -I _i_ lit — IL IL — III -~ , IIL-f- IL 

B p 2 [xq~^~p~^~ ) B p2 {x~ 1 p ■ Li ~ q-^i— ) fc) p2 [xq~^~p 



m-\-n m-|-n > 

2 



(5.13) 



m-f-Ti -1 
ra+m m-|-ra » 9 

2 



p2 (xp 1 2 g 

/ 1 1 ra + m m + n . jib 

6 p 2(x-V 2 g 2 ) i=i+rL 



n 



6^2 (x Vp O^iM^P 1 V) 



2 P 



Q^NV 1 ) ^ 1 *p x V) 



Albeit this formula has been obtained for both n and m even, it is easy but lengthy, by 
following the same path, to show that indeed formula ()5.13|) gives the scattering factor 
between the general n-th and the m-th breather. However, two precautions should be 
observed when dealing with (|5.13|) . The first one is that the product over / is present 
only if m > 2; the second one concerns that this product runs over the integers if n 
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and to have the same parity, whereas over the half-integers if n and to have different 
parities. We end this part by noticing two simple properties of (|5.13|) 



Sn,m(@)S nt m( 9) — 1 , S n>m (9) — S m;n (9) . 



(5.14) 



Remark 1: The scattering factors between the elliptic breathers, ()5.13j) . are not 
independent. More precisely, they satisfy the decomposition 



S, 



n+m,n'+m' 







iK iK< 



(n — to — n' + to') 



Sm,m' 
Sn,n' 



o- 


(iK 




e- 


(iK 




H 








(iK 





2 2r] 

[n + to + n' + to') 



2r/ 
iK' 
~2rj 
iK 2 
~2r[ 
iK 2 
~2rj~ 



(— n — m + n! + m!) 
(n + to — n' — to/) 
(— n — m — n' — to') 



(5.15) 



In particular, this links the general scattering amplitude S^n to that of the fundamental 
scalar particle Sii (or, in mathematical language, to the DVA structure function): 



n-l 



S n , n (9) = Y[ S hl n -W 



l=l-n 







V 



(5.16) 



Remark 2: Let B n (9) be the Zamolodchikov-Faddeev generator of the n-th breather 
with rapidity 9 and satisfy the exchange algebra 



B n {9\)B, 



m\V2, 



s 



9o)B 



n, 



B„ 



(5.17) 



The relation ()5.15|) is consistent with the following connection between the creation 
operator B n+m {9) and the operators B n {9) and B m {9): 



B„ 
B,, 



' (iK iK 2 \ . 
9 + (to - n) 

(n + TO) 



„ , «K *K 2 . 



(5.18) 



Consecutive applications of this property allows us to write the generic breather gener- 
ator B n as a product of fundamental scalar generators B\. This fact becomes important 
as soon as we realise [5] that B\{9) coincide with the mode generator of DVA, T(z), 
which enjoys a particular structure pp. Hopefully, this should clarify the physical role 
played by DVA. Finalising the parallelism, (|5.18|) is the elliptic deformation of the 
analogous relation in the sine-Gordon theory (see footnote 9 in the seminal paper [H]). 
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5.1 Recovering the trigonometric limit 

Let the trigonometric limit be defined by K' — ► oo or, equivalently, by K — > tt/2. In 
this respect a key limit is 



Q p 2(x qp 

e p2 (x- i gy- z ) 



1 



tanh 



(5.19) 



where we have introduced the quantity b defined by r], and then, thanks to the nature 
of that limit, connected to the sine-Gordon coupling constant (see section 2) as follows: 



7T 



rj 



ff 1 



(5.20) 



2(6 + 1)' l-/? 2 ' 
Application of the limit ()5.19j) to the scattering factor (|5.1Hjl gives for us 

1 tanh 2 (f + fZ) 



tanh [| + if (n - m)] tanh [§ + if (m + n 
tanh [§ - if (n - m)] tanh [f - *f (m + n 



2 



n 



tanh 2 (| - if Z) 



sinh # + z sin ^{n — m) sinh 6 1 + z sin f (m + n) 



sinh # — i sin ^{n — m) sinh — i sin ^(m + n 



-b i 



m—l 



n 



sin 2 [if - f (n - m + 2/)] cos 2 [if + f (n - m + 2Z)] 



J sin 2 [-if - f (n - m+ 2/)] cos 2 [-if + f (n - m + 21)] ' 



(5-21) 



As expected, this limiting formula coincides with the scattering amplitude between the 
n-th and the m-th breather in the sine-Gordon model (i.e. formula (4.19) of [0]). 



6 Higher Order Deformed Virasoro Algebras 

Continuing the identification between the fundamental scalar generator Bi(8) and the 
DVA generator T(z), we can now define higher order Deformed Virasoro Algebras. 
The order n deformed Virasoro algebra is the central extension of the Zamolodchikov- 
Faddeev algebra 2 of two breathers creators B n , i.e. 



B n (6i)B ri 



S 



n,ri\ u l 



0o)B r , 



)B 



(6.1) 



where 



S n ,n{@) — x 



O p 2(x 1 q n p 



)Q p 2(xp 1 - n q n ] 



Q„2 (xq n p~ 



)Q p 2(x- 1 p 1 - n q n ] 



n-l 

n 

i=i 



X 



,0 2 2 (x 1 q l p l )QK(xp 



%{xq l p- l )Q 2 ^{x-^- l q l ) 



(6.2) 



When n = 1, the algebra ()6.1|) coincides with the Deformed Virasoro Algebra Vir Piq 
jTj without central extension. In order to clarify the relation between the general ()6.1|) 

2 Some authors prefer to define the Zamolodchikov-Faddeev algebra with the central extension: in 
this case it would simply coincide with a higher order deformed Virasoro algebra. 
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(or, which is the same, the corresponding higher order deformed Virasoro algebra) and 
the particular Vir PA , we consider the semi classical limit of ()6.1|) around the particular 
values q = p h (or K = rjh), h = 1,2, characterised by the special property that at 
them S n ^ n (6) = 1. From the algebraic perspective, this means that we can define a 
Poisson bracket around each point for which q = p h . In this respect, the order n 
Zamolodchikov-Faddeev algebra (generated by B n ) turns out to be a quantisation of 
such Poisson structure. Explicitly, if we measure the distance through f3 such that 



p h , these Poisson brackets are defined by 



{B n {6 x ) , B n {6 2 )} 



d_ 

dp 



S. 



n,n\ u l 



02)\f3=oB n (9l)B Tl 



where 







8(3 

oo 



Sn,n(0l2)\/3=Q = C(n, fl) hip 



X 



X 



1 — X 



-1 



m=0 



2x p 



1^2m 



x -lp2m 



+ 



2x p 



l„2m+l 



1 — X 

,2m 



x -l p 2m+l 



+ 



2xp 



2xp 



2m+l 



1 — xp 



2m 



1 - xp 2m+1 



(6.3) 



(6.4) 



and besides 



c{ n ,h)= 



(l-h)l 



(n 



l=l-n 



h 

'2' 



\h+l 



C(n,l) 



n 
~2 



C(n,2) 



'-IY-1 



(6.5) 

Interpreting this relation in terms of modes (see for instance the discussion in Section 
5 of ^H])> we remark that the Poisson structure around q = p h is proportional to the 
Frenkel-Reshetikhin Poisson bracket [TOj, which defines the classical Deformed Virasoro 
Algebra, without central extension. Even this conclusion has led us to consider the 
breather algebras as higher order versions of the fundamental DVA 3 . 



7 Conclusive outlook 

In this paper we have continued the study of scattering factors between bound states 
of the spin 1/2- XYZ model, commenced in [H] and restricted there to the simplest one. 
Following the method exposed in [H], any scattering factor has been extracted from the 
infrared limit of the NLIE describing the pairwise collision. In the end, we have found 
the complete set of scattering factors between bound states (elliptic breathers) of the 
spin 1/2-XYZ chain (in the attractive regime). This physical information has allowed 
us to re-interpret the braiding relations between the corresponding scattering operators 
as a consistent mathematical definition of a family of Deformed Virasoro Algebras, 

3 Notice that it is 'fundamental' even in a physical sense as it describes the scattering of the funda- 
mental scalar particle. 
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parameterized by a positive integer n. Of course, at n = 1 the usual Deformed Virasoro 
Algebra [1] is recovered and n ranges from this value to a maximum depending on the 
chain parameters. And all the members of this family have been proved to reproduce, in 
the classical limit, the Poisson bracket structure discovered by Frenkel and Reshetikhin 
[TO] (and defining the classical Deformed Virasoro Algebra). For these reasons, these 
algebras are all eligible to the name of higher order Deformed Virasoro Algebras and 
generalise the seminal quantisation by Shiraishi-Kubo-Awata-Odake pp. Besides, the 
search for a family of DVAs seems to be a quite intriguing problem, though we cannot 
guaranty that we have found the amplest family neither the only one. Nevertheless, 
these results should shed light on the peculiar role of DVA in field and lattice theory. 

In an incoming paper [T2] these results and others will be used to highlight the 
non-relativistic nature of the thermodynamic theory and the link between bound state 
energies and the poles of the Zamolodchikov Z4 S-matrix 
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